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Abstract. In this paper, we apply the thermodynamic framework recently put into place by 
Rajagopal and co-workers, to develop rate-type models for viscoelastic fluids which do not possess 
instantaneous elasticity. To illustrate the capabilities of such models we make a specific choice 
for the specific Helmholtz potential and the rate of dissipation and consider the creep and stress 
relaxation response associated with the model. Given specific forms for the Helmholtz potential 
and the rate of dissipation, the rate of dissipation is maximized with the constraint that the 
,_i difference between the stress power and the rate of change of Helmholtz potential is equal to the 

rate of dissipation and any other constraint that may be applicable such as incompressibility. 
We show that the model that is developed exhibits fluid-like characteristics and is incapable of 
instantaneous elastic response. It also includes Maxwell-like and Kelvin- Voigt-like viscoelastic 
materials (when certain material moduli take special values). 
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1. Introduction 



Recently, a thermodynamic framework has been put into place to describe the response of 

dissipative bodies that includes a large class of viscoelastic bodies (see [1] and [2], for details of 

>. the framework). With regard to the response of viscoelastic bodies, they consider the response 

to be that of a class of elastic bodies from an evolving set of configurations which they refer to as 

natural configurations (also see [31 II])- The evolution of the natural configuration is determined 

C<~) by the rate of dissipation, or to be more precise, the maximization of the rate of dissipation. 

r~-- In a purely mechanical context, the response of the material is characterized by constitutively 
prescribing the stored energy (or Helmholtz potential) and the rate of dissipation functions. 

,_i Since in a closed system the entropy increases to achieve its maximum equilibrium value, the 
quickest way in which the maximum could be reached is by maximizing the rate of dissipation. 
Of course, this is a plausible assumption and not a "principle" . Also, to require such a criterion 
for an open system is not on very sound footing. However, it is surprising how well such a 
c3 requirement works. Using such a thermodynamic framework a variety of material responses such 
as viscoelasticity [HI E], twinning [TJ, solid-solid phase transition [8], plasticity [9], crystallization 
of polymers [TUl ITT] , single crystal super alloys [I21IT3"] . response of multi-network polymers [H] 
and anisotropic liquids [T5] have been modelled. Particularly in viscoelasticity, this framework 
has been used to generalize one dimensional models due to Maxwell [16], Kelvin [17] and Voigt 
[IS] . Burgers [19], and the standard linear solid to three dimensions [20l EH G2]. Moreover, 
recently, it has been shown within the context of Maxwell fluid [23] that by choosing different 
forms for the stored energy and the rate of dissipation, one can obtain the same constitutive 
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relation for stress. It has also been shown (see [24j) that this framework leads to more than 
one three dimensional model that reduce within the context of one dimension to Burgers' fluid 
model. 

Based on the aforementioned general framework, Rajagopal and Srinivasa [5] have developed 
a methodology to obtain rate-type models for fluids which possess instantaneous elasticity and 
in particular have derived a three dimensional generalization for the one dimensional Maxwell 
model (with the mechanical analog - a spring and a dashpot in series). When one assumes 
that the displacement gradients are small, this non-linear model leads to the classical upper 
convected Maxwell model. The aim of this paper is to use the framework to develop rate-type 
models for fluids which do not possess instantaneous elasticity. Specifically, we shall assume 
that the response from the natural configuration to the current configuration is akin to that 
of a generalized Kelvin- Voigt solid. Upon removal of external load, the body moves to the 
natural configuration with some "relaxation time" which is greater than the intrinsic time t m 
(also see [25]). By setting this relaxation time to a value less than t m , one can obtain the class 
of models that can be generated using the framework in [5] and thus our work in this paper can 
be viewed as a generalization of the analysis in [5]. As an example, using our framework, we 
develop constitutive relations which in one dimension reduces to a dashpot and a Kelvin- Voigt 
element (which is a spring and a dashpot in parallel), in series (see Section [5]). In carrying out 
the maximization of the rate of entropy production, one needs to decide what one maximizes 
with respect to. In most of the studies that have been cited, the maximization is with respect to 
appropriate kinematical variables and determining a constitututive representation for the stress. 
On the other hand, one could maximize with respect to the stress. In all the problems considered 
thus far, both methods lead to the same answer in explicit constitutive theories. However, it 
is possible to obtain a much larger class of constitutive relations following the latter procedure 
(see [26]). Also, from a philosophical perspective it is preferable to use the latter procedure. We 
shall not get into these issues here. In this paper we maximize with respect to the kinematical 
variables and obtain a constitutive expression for the stress. 

While the idea of maximizing the rate of dissipation was also enunciated by Ziegler [27] , Wehrli 
and Ziegler [2H] our philosophy, interpretation and use of the maximization of the rate of entropy 
production is quite different, and the difference is not minor. We have discussed the differences 
in our approach in some detail in [26]. In fact, Ziegler's approach cannot be used to obtain a 
whole host of models including a whole class of implicit constitutive relations (see [26]). Also, 
the original work of Ziegler contains certain mathematical errors as pointed out in [2] . 

After a discussion of the preliminaries in the next section, we develop a general framework for 
the development of constitutive models for viscoelastic bodies that do not possess instantaneous 



elastic response in Section 3.1 Within our general framework, we derive a specific model in 



Section 3.2 which stores energy like a neo-Hookean solid with a rate of dissipation which depends 



on the stretching tensor of the natural configuration and the stretching tensor between the natural 



configuration and the current configuration. In Section 3.3, we show that our model reduces 



to either the Maxwell-like fluid or the Kelvin- Voigt-like solid under certain restrictions on the 



material parameters. We shall solve the problem of uniaxial extension in Section |4.1[ followed 
by some remarks on the application of our framework to develop models for visco-elasto-plastic 
response in Section [6] 



2. Preliminaries 

Let Kr(B) and n t (B) denote the reference configuration and the configuration of the body B 
at time t (or the current configuration), respectively. Let X denote a typical point belonging to 
k r {B) and x the same material point at time t, belonging to n t {B). Let Xk r denote a one to one 
mapping that assigns to each X G kr{B), a point x G Kt(B), i.e., 

x:= XKfl (X,t). (2.1) 

We shall assume that Xk r is a sufficiently smooth mapping. The velocity v, the velocity gradient 
L and the deformation gradient F Kfl are defined through 



m ' ax' K * ' ax 



v:=^, L:=^-, F Kfl :=^. (2.2) 

It immediately follows that 

L = P« H F~i, (2.3) 

and the symmetric part of the velocity gradient D is given by 

D := 1 (L + L T ) , (2.4) 

where (.) T denotes transpose of a second order tensor. The left and right Cauchy-Green stretch 
tensors B Kfl and C KR are defined through 

B KR := F KR F KR , C KR := Y Kr F Kr . (2.5) 

Let K p (t) denote the natural configuration associated with the configuration n t . We define F K 
as the mapping from the tangent space at a material point in k p u) to the tangent space at the 
same material point at n t (see Fig. [I]). Similar to Eq. (2.5), we can also define the left and right 
Cauchy-Green tensors from the natural configuration to the current configuration n 

T3 Tp TpT /-( TpT Tp 

""■Pit) -~ r K p(t) r K p(t) ' ^ K p(t) " Kp{t) K PW 

The mapping G is defined through (see Fig. [I]) 

C •— P 1 ■= P 1- -^ P" 

^ • r BH->K p (,) • «p(t) K R' 

We define the tensors B G , L G and D G through 

B G := B Kfl _s. Kp(t) = GG , L G := GG , D ( 

In addition, let the tensors L p ( t ) and D p ( t ) be defined by 

L p (t) := F Kj(t) F^ w , D p(t) := - (L p(t) + Lj (t) ) . (2.9) 

Also, the principal invariants of B p /^ are denoted by 

Ib pW = tr(B p{t) ), II Bp{t) = i [(tr(B p(t) )) 2 - tr(Bj (4) )] , J// Bp(t) = det(B p(t) ), (2.10) 

where tr(.) is the trace operator for a second order tensor and det(.) is the determinant of a 
second order tensor. 

In this paper, henceforth we shall suppress k and denote B K by B p ( t \. 
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(2.6) 




(2.7) 


5 (Lo + Lg) . 


(2.8) 



Now, from Eq. (2.7): 



r Kn — r K ... Cjt H - F K ... tc 

K -R re J>(*) K p(*) 



F F _1 



F Kp(t) GG 1 F Kp 1 (i) +F Kp(t) G 



(2.11) 



where 



^L-L pit) + F Kp(t) L G F K ^ t) . 

is the material time derivative of the second order tensor given by 

0(A) 



A: 



Of 



+ grad(A)[v] 



(2.12) 



for a second order tensor A with grad(.) being the gradient of a second order tensor with respect 

(2.13) 



to the current configuration K t . Hence, from Eq. (2.11), 

tT 



p(t) "r r « p (0 G *p(*) ' 



and so 



D = D 



p(«) + 2 



f^LgF: 1 +f: t .l^f t 



S(t) K p(*) 1J ^ !J 



Pin 



(2.14) 



Further, we note that the upper convected Oldroyd derivative [29] (Also see [9] for the inter- 
pretation of Oldroyd derivative within the context of theory of multiple natural configurations) 
of B p (i) can be related to Dq through (see [5]) 



v 

B P (t)--- 



B 



Pit) 



LB 



P(t) 



B p (t)L 



-2F K ,..~D n F T 



K p(t) 



'G* 



>(t) 



(2.15) 



Assuming that the body under consideration is incompressible, we shall record the balance of 
mass, and the balance of linear and angular momentum (in the absence of body couples): 



div(v) =0, pv = div(T T ) + ph, 



(2.16) 



where p is the density, v is the velocity, T is the Cauchy stress tensor, b is the specific body 
force and div(.) is the divergence operator with respect to current configuration n t . In addition, 
the local form of balance of energy is 



pe = T.L — div(q) + pr, 



(2.17) 



where e denotes the specific internal energy, q denotes the heat flux vector and r denotes the 
specific radiant heating. 

Finally, we shall assume that the body under consideration undergoes isothermal processes. It 
is easy to modify the procedure to take into account non-isothermal processes. We shall invoke 
the second law of thermodynamics in the form of the reduced energy dissipation equation, given 
by (see [30112]): 

T.D-p^ = p#C:=£>0, (2.18) 

where T is the Cauchy stress, ip is the specific Helmholtz free energy, ( is the rate of entropy 
production and £ is the rate of dissipation]^ 



The terminology rate of dissipation usually refers to the conversion of working into energy in thermal form 
(heat). However, while considering general non-isothermal processes one uses the term to the product of the 
density, temperature and entropy production. 



3. Constitutive assumptions 

3.1. General framework. In this subsection, we shall first constitutively specify general forms 
for the specific Helmholtz potential ip and the rate of dissipation £. Using Eq. (2.18) and incom- 
pressbility as constraints, we shall maximize the rate of dissipation £, to obtain our constitutive 
relations. We shall assume that during its motion from K p (t) to Kt, the body stores energy as well 
as dissipates. The storage is due to elongation of the polymer networks. Now, assuming that 
the response of the body from K p (t) to K t is that of an isotropic body along with the assumption 
of incompressibility, we shall choose the specific Helmholtz free energy to be a function of the 
Cauchy-Green left stretch tensor B p ( t ), i.e., 



This implies that, 



v> 



i> = ^( b p(*)) = 

dip dip 



ip{lB Kt) JlB pit) ] 



dh 



+ h 



'p(t) 



dlh 



B 



p(t) 



dip 



dlh 



T)2 



D 



P (ty 



(3.1) 



(3.2) 



Also, we shall assume the rate of dissipation to be a function of the stretching tensor and Cauchy- 
Green left stretch tensor between the K p i t ) to K t , the stretching tensor between k r and K p ( t ), i.e., 

£m = (m(B p (t), D p ( t ), D G ). (3.3) 

In other words, the body dissipates energy during its motion from K p (t) to Kt (due to continuous 
scission and healing of polymer networks) and also dissipates energy during its motion from Kr 
to K p (t) (due to sliding of polymer strands over one another). 

Due to the assumption of isotropic elastic response, the stored energy remains unchanged 
under any rotation. Hence, for our calculations, we shall assume that the natural configuration 
is rotated such that U 

K p(t) = ^ K P (t)' (3.4) 

and therefore 



D 



p(t) 



V V -1 +V _1 V 



rV: 1 b pW v 



>(*) 



>(*) 



On substituting Eq. (3.2) into Eq. (2.18), we get 

T.D — Tp( 4 ).D p ( t ) = £ m ( 

where 

dip dip 



Bp(t),D p ( t ),D( 




dip 



dll 



B 



- B l(t) 



Substituting Eq. (|2.14|) into Eq. (|3.6|), we obtain that 

2 T - [ Fk pW^- 
Since T is symmetric, Eq. (3.8) reduces to 

F^TF K -J t) .L G +(T-T p(t) ).D p(t) =^ 



p(t) 



vBp(t),D p ( t ),D 



G) 



(Bp(t),D p ( t) ,D 



G) 



(3.5) 

(3.6) 
(3.7) 

(3.8) 
(3.9) 



For the application of this thermodynamic framework to anisotropic fluids the reader should refer to [15] 
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Also the assumption of incompressibility leads to 

tr(D) = tr(D p(t) ) = £r(D G ) 



0. 



(3.10) 



Now, following Rajagopal and Srinivasa [5j, we maximize the rate of dissipation £ m by varying 
L G , Dp( t ) for fixed B„( t ) with Eqs. (3.9), (3.10) as constraints. We maximize the auxiliary function 



~Zrn-Fl (t TF^ t) .L G -(T-T p{t) ).B 



$ given by 

$ := U + Ai 

+ A 2 (I.D G ) + A 3 (I.D p(t) ), 
where Ai,A2,A3 are Lagrange multipliers. By setting, d&/d~D p ( t ) 



T p (t) + ^-1 



Ai 



d~D 



p(t) 



and 



A, 



I 



Ai + 1 



,— T ^« 



P(4) J (3.11) 

and d&/dL G = 0, we get 

(3.12) 

(3.13) 



Ax ' V Ai J Kp «<9L G K *w 
At this juncture it is worth recalling the comments in the introduction concerning maximization 
with respect to the kinematical quantities. On substituting Eqs. (3.12), (3.13) in Eq. (3.9), we 
finally obtain 

'Ai + 1\ £„ 



Ai 



ggg 



d£r. 



9L G - L G + OD 



D 



(3.14) 



Pin 



p(t) 



Hence, from Eqs. (|3.12|), (|3.13|), the final constitutive equations reduce to 

dip 



pl + 2p 



01 



+ h 



B 



pit) 



dip 
' p(t) dII B . 



B 



dip 



v(t) 



pW 



dll 



B 



B 



Pit) 



pit) 



+ 



Zr, 



di, 



9£,m T I 



d£,„ 



9D 



P(t) 



D 



p(t) 



<9D p(t) ' 



(3.15) 



AI + 



ae 



<9L G ' 



• L G + 



£>£„ 



9D 



p(t) 



D 



p(t) 



sw 9L G 



m -rpT 



>(*) ' 



where A := y-,p '■= t 1 are Lagrange multipliers due to the constraint of incompressibility. At this 
point it appears that there are two constitutive relations - (3.15 )„, (3.15 )& - for stress instead of 
just one! We would like to note that the two expressions - (3.15 ) , (3.15 )& - have to be equated 
and simplified to obtain one expression for stress and an evolution equation for B p ( t ). This will 
become clear within the context of the specific choices of ip and £ that are made in the next 
subsection. Furthermore, expression (3.15 )& appears at first glance non-symmetric but is in fact 
symmetric for the specific form chosen for £ as we shall show in the next subsection. 

3.2. Specific model. We now derive constitutive expression for the stress by choosing the stored 
energy to be 



Hh 



p(t)i B p(t)< 



2p 



/b pM " 3 



(3.16) 



and rate of dissipation to be of the form 

£m(Bp( t ), D p(t ), D G ) = ?7pDp( t ).B p ( t )Dp( t ) + ?? G D G .Bp( t )D G . (3.17) 

The stored energy chosen here is that for a neo-Hookean material with \x being its elastic modulus, 
whereas the rate of dissipation is similar to that of a "mixture" of two Newtonian-like fluids (in 
the sense that the dissipation is quadratic in the symmetric part of velocity gradient), whose 
dissipation also depends on the stretch (specifically the stretch from the natural configuration to 
the current configuration), with viscosities r/ G and ^Jj The former term on the right hand side of 
expression (3.17) is due to the dissipation during the motion from k r to K P (t) and the latter term 
is due to dissipation during the motion from k p u\ to K t . Note that with the above choices for ip 
and £, as the body moves from the k p m to « t , there is both storage (like a neo-Hookean solid) 
and dissipation (like a Newtonian-like fluid) of energy simultaneously and hence, k p m evolves 
like the natural configuration of a Kelvin- Voigt-like solid (also see [20]) with respect to K t . Now, 
with this choice for ijj and £, the constitutive relations given by Eq. ( |3.15 ) reduce to 

(3.18) 
(3.19) 



77 

pi + /xBp( t ) + -^ (Bp( f )D p ( t ) + Dp( f )Bp( f) ) , 



and 



AI + f P^ t) (B p(i) D G + D G Bp (t) ) F T Kp(t 



Now, from Eq. (3.18) and Eq. (3.19) 



\p - A J I + /iBp( t) +y (Bp( t )Dp( t ) + Dp(i)Bp( t) ) 



2 k pW 



(B p (t)D G + D G B p(t) ) V Kp(t) , 



(3.20) 



and using Eq. (2.15) and Eq. (3.5) in Eq. (3.20), we get 



(p - a) i+^ B p (t) + 1 (v Kp(t)B p (i) v K -; (t) + v;; (t) Bp (t) v Kp(t 

v _ ± _ 1 v 

K P(t) Bp(t) V k pW "" ^K p(t) Bp(i) V Kp(t 



Vg 



-V" 1 
4 K p(t) 



V 



K p(t) ' 



(3.21) 



Pre-multiplying and post-multiplying Eq. (3.21) by V K , we have 



P 



> P (i)+A i B 

no 



V P 



;■<*) + 4 



Bp(t)Bp (t) 



Bp(t)Bp (t ) 



V" 1 

I k p(*) 



v v 

-Bp( t ) Bp(t) - Bp(i) Bp( t ) ) V Kp(i) . 



Also, fromEq. (3.21) 



V 



A 



yfr(Brf.)J+^tr( B /<<'> 



/^ r ( b p(<)) 



Eq. (3.19) can be re- written as 

T = AI + —V" 1 

4 K vw 



-l -l 

~ V «p( t ) Bp(<) V K p(t ) ~~ V « P( t) B p(*) Vk pm 



V, 



pW 



(3.22) 



(3.23) 



(3.24) 



Of course, one can choose the rate of dissipation to be quadratic in the symmetric part of velocity gradient 
without any stretch dependence, for example, £ = rjpDpM.TipM + 7/gDg.Dg. The resulting model would be a 
variation of the current model. See Section [71 for the calculations using such a rate of dissipation. 
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Eq. (3.22) and Eq. (3.24) are the final constitutive relations with Eq. (3.22) togther with 



Eq. (3.23) being the evolution equation for the natural configuration K p (t). 

In what follows, we shall also show that our model can be reduced to both the Maxwell- 
like fluid model and Kelvin- Voigt-like solid model under certain assumptions for the material 
parameters that are involved. We shall solve the problems of creep and stress relaxation by 
considering homogeneous uniaxial extension. Based on the results for creep and stress relaxation, 
and following the definitions given in J3TJ [32] for a fluid-like body and a solid-like body, we shall 
show that our model is a fluid-like model when none of the material parameters are ignored. 



3.3. Limiting cases. By setting r/ p to zero, one can see from Eq. (3.20) that B p w an d Dg have 
the same eigen- vectors and hence commute. Using this fact, Eq. ( 3.18p and Eq. ( 3.22[ ) reduce to 
the Maxwell-like fluid model developed by Rajagopal and Srinivasa (see Eqs. 40-42 in |5J). On 
the other hand, if we assume t\q goes to infinity, with i] p , fi and all other kinematical quantities 
remaining finite, then as the deviatoric part of Eq. (3.19) is finite, from Eqs. (3.19), (3.10), we 



must have D G — > 0. This also implies that G is pure rotation and hence B p ( t ) -> B Kfl , D„( t ) — > D. 



Thus, Eq. (3.18) reduces to 



T = pi + ^B K 
This is a generalized Kelvin- Voigt solid model. 



+ | (B KB D 



DB„ 



(3.25) 



4. Application of the model 

4.1. Creep. In this subsection we shall solve the problem of homogeneous extension under con- 
stant stress and show that our model is a fluid-like model when none of the material mod- 
uli are ignored. Before we go into the details of the problem, we would like to note that for 
steady problems wherein the deformation is homogeneous, the results for our model would be 
same as that for a Maxwell-like fluid with stretch dependent relaxation as formulated in [5]. 
This is because B p ( t ) = (since —^p- = 0, due to the assumption that the deformation is 
steady and grad(B p ^) = as the deformation is homogeneous) and the constitutive relations in 



Eqs. (J3T22J), fl3T24[ ) reduce to Eqs. (36), (39) in [5]. 

Now, in the case of time dependent homogeneous extension: 

x = X(t)X, y 



-Y, 



the deformation gradient with respect to the reference configuration is given by 



K R 



diag< \{t), 



VWiWWJJ 



Hence, the velocity gradient is given by 



diag 



A A 



A 



A' 2A' 2A 



(4.1) 



(4.2) 



(4.3) 



Now, we shall assume that 



B 



pit) 



diag^(t), 



y/wy vw) 



(4.4) 



and hence 



B P (t)= diag<^ B(t) 



B p (t) 
2B(t)X(t) 



, , , Bit) Bit) } 



B(t) 



+ 



X(t) 



B(t) 



X(t) ' 2BW(t) VB(t)X(tY 2BW(t) 



+ 



_A(*) 1 
B(t)X(t)j 



and 



V, 



p(<) 






diag i VB(t) '^W'^W/ 



(4.5) 
(4.6) 

(4.7) 



For the case of homogeneous extension, Eq. (3.18) and Eq. (3.19) reduce to 



T 
T 



pi + /xBp(i) + T)pB p (t)Dp(t), 
XI + r/ G Bp( t )D G , 



and the final constitutive relations Eq. (3.22) and Eq. (3.24) become 



-i 



AI - ^-V Kp(t) B pW v- j 



(*)' 



(p - A)I + /xBp( t) = — — B p (t) - yB p ( t) , 



where 



(p-A) 



3// 



tr(B 



p(t)' 



On substituting Eqs. (4.5), (4.6) into Eq. (4.10), we arrive at 

-B 5 / 2 (t) - Bit) 



Bit) 



2 VG B(t)X(t) 4/i 



Vg + Vp A(t) 



Vg + Vp 



l + 2B 3 / 2 (t) 



(4.8) 

(4.9) 
(4.10) 

(4.11) 
(4.12) 



Using Eqs. (4.6), (4.7) in Eq. (4.9) and using the fact that lateral surfaces are traction free, we 
find that 

1 \ (2B(t)X{t) 



T u 



2 



1 + 



Bit) 



2B 3 / 2 {t)J \ X(t) 

Solving Eqs. (4.12), (4.13) simultaneously, we obtain 

2Tn 4/i / B 5 / 2 (t) - B(t) 

^ V l + 25 3 / 2 (t) 



B(t) = - 

^P [} + 2B3/2(t) / 

A(£) 2(77p + 77 G ) / T lt B l / 2 (t) 



x(t) ~ 

which can be re-written as 



^p 



2ju / fi 3 / 2 (t) - 1 \ 



dB 



1 + 



l + 2B s / 2 (t)J r] p \l + 2B 3 / 2 {t)J 
2T n _jB*l\t)-B{ty 



2S3/ 2 (t) 



1 + 25 3 / 2 (i) 



IdX 

x ~M 



n (l Xf fnB^jt) \ ( B*/ 2 (t)-l \ 

\rj J \l + 2B 3 / 2 (t)J \l + 2B 3 / 2 (t)J 



(4.13) 

(4.14) 
(4.15) 

(4.16) 
(4.17) 



where t 



f u 



- Jx is a non-dimensional stress, 



and 77 



— is a non-dimensional time 
ratio of viscosities (r]p,r] G ). With known values for T n and fj, the ODEs - Eqs. (|4.16|), (|4.17|) 



^ is the 

VG 



were solved simultaneously using "ode45" solver in Matlab. The initial conditions B(Q) = A(0) = 
1 were used for the loading process. In the case of unloading, the initial values were set to the 
values of B, A evaluated at the end of loading process. Before we discuss our results, the reader 
should note that B is the square of the stretch from the natural configuration to the current 
configuration and A is the stretch from the reference configuration to the current configuration 
or the total stretch. 

We shall now discuss the numerical results obtained for the problem of creep. Fig. pi) portrays 
the results (specifically, B as a function of i and A as function of t) for Tn = 1 for various values 
of fj. For all the cases of fj, the evolution of the natural configuration with respect to the current 



configuration (or Basa function of t) is the same. This can also be seen from Eqs. (4.16), (4. 17) 



as the evolution of B in time depends only on the value of Tn, and does not depend on fj. Also, 
as seen in Fig. pi) (a), on removal of the load, the stretch in the body does not return to unity 
but exhibits a permanent residual stretch, which is the typical behavior of a fluid-like body (see 
[311 [32]). I n addition, from Fig. (J2J) and Fig. (J3J), by increasing Tn (from 1 to 5) for fixed fj (here 
the value is 10) - the maximum value of B increases, the maximum value for A increases and the 
permanent residual stretch also increases. 



Now, as discussed earlier in section 3.2, the natural configuration evolves like a Kelvin- Voigt- 
like solid with respect to the current configuration. Fig. pi) (b) reiterates this fact since B varies 
with t in a fashion similar to the stretch as a function of time for a Kelvin- Voigt solid in a creep 
experiment. Fig. pi) (a) also shows that, as fj increases, the maximum value for total stretch (A) 
decreases. This shows that for a fixed amount of loading time (for instance for t = 10) as the 
value for fj increases, the value for maximum total stretch decreases and so the rate of relaxation 
decreases. As fj increases one would expect that our model would tend to a solid model in the 
limit of fj —?• 00. To see this, in our numerical simulations, we set fj to a very large number 
(specifically f\ = 100000), with T n = 1, and we found yB and A to follow the same trend in 
time (see Fig. 111). This is an extreme case when the natural configuration and the reference 
configuration tend to being the same, reducing our model to Kelvin- Voigt-like solid model. 

4.2. Stress under constant strain rate. If we define e := ln\ as our strain measure, then 

■ B 5 / 2 (t) - B(ty 



Eq. (4.12) reduces to 



B{t) 



a 



_2f]__ 
Vg + Vp 



B(t)e 



4yU 

VG + Vp 



Using Eq. (|4.18|) in Eq. d4.13|), we get 

1 



T 11 = r iG [l + 



2B*l\t)J Iv p + Vg 



'Op 



-B(t)e 



1 + 25 3 / 2 (t) 

2p 5 5 / 2 (£) - B{t) 
r] G + r] p 1 + 253/2 (t) 



Upon non-dimensionalizing Eqs. (4.18), (4.19), we arrive at 



dB 2fjB de 4 / 5 5 / 2 - B \ 

~dl ~ fj+l'dl ~ 77 + I \1 + 2B 3 / 2 ) 



Tn 



253/2 



Brj de 



( 



55/2 _ B 



1 + fjdt I+77 V 1 + 253/2 
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(4.18) 

(4.19) 

(4.20) 
(4.21) 



With known values for 4=, B can be evaluated using Eq. (4.20 ) and then from Eq. (4.21 ), T n can 



be calculated. For stress relaxation, we set % to zero and solved the ODEs - Eqs. (4.20), (4.21) 
- numerically using the solver "ode45" in Matlab. Fig. (pi) displays the plots for B, and T n as 
functions of i for various values of fj. The initial condition for B was chosen to be 2. Fig. (|5])(b) 
also shows that the stress finally relaxes to zero. This is a characteristic of a fluid-like material. 
In addition, one can also see from Fig. d5J) that as the ratio of viscosities fj increases, the time 
of relaxation increases. Also from Eq. (4.21), as fj — > oo, Tn — > and hence cannot stress relax. 



This is because as i] — > oo, the model behaves like a Kelvin- Voigt-like solid model as can be 



gleaned from the creep response discussed in Section (4.1) . In addition, one can also notice from 
Fig. (|5])(b) that the initial value for T n decreases as fj increases. 



5. Model reduction in one dimension 
If the elastic strain is small in the sense that 

||B p(t) - I|| = 0(e), e<l, 



then Eq. (4.11) reduces to (see [5]) 



P- A =-!*, 



(5.1) 
(5.2) 



and hence Eq. (4.8) reduce to 



(5.3) 



T = AI + r/ G B p(t) D G , 

H (B p{t ) ~ I) = B p(t) (r/ G D G - ^D p(f) ) . 

If Aj (i = G,p) is the one-dimensional stretch and e = InXi is the true strain, then in one 
dimension, Eq. (5.3) reduces to 

. A G 



a 



Vg 



n(x 



Xg' 
1) 



Ar 



A, 



(5.4) 



t)gt- - %t-, 

A G A p 



where fji = r/jA? (i = G, p) are the stretch dependent viscosities and a is the one dimensional 



stress. Eq. (5.4) under the assumption that e^ <^i 1 

2fi€ p = 



a 



Vc^g, 



(i = G,p) reduces to 
Vg^g - Vptp- 



(5.5) 



Eqn. (5.5) can also be obtained if we have a dashpot with r} G as the viscosity and a Kelvin- Voigt 
element, with a spring constant 2/x and viscosity fj p , in series (see Fig. fl|. Thus, our model 
in Section |3.2| reduces to a dashpot and a Kelvin- Voigt element in series whose viscosities are 



stretch dependent. 



6. Concluding Remarks 



For the model developed in this paper, we solved the problems of creep and stress relaxation 
and based on the results showed that our model is a fluid-like model. We have also shown that 
under certain idealizations, our model reduces to both a Maxwell-like fluid and a Kelvin- Voigt- 
like solid. We have also shown that our model is a three dimesional generalization of the one 
dimensional model with a dashpot and Kelvin- Voigt element in series. We would also like to 

n 



note that the general framework developed in Section 3^ can be extended to model visco-elasto- 
plastic response by choosing the rate of dissipation for the motion from k r to K p (t) to be of the 
form given in Eq. (39) or Eq. (40) in [5]. For example, one can choose the rate of dissipation to 
be of the form 

1 /9 

fm(B p (t), Dp( t) , D G ) = Y (D p(t) .Bp (i) D p(t) ) + 77 G D G .B p(4) D G , (6.1) 



where Y is a material constant, along with the stored energy in Eq. (3.16), and then maximize 
the rate of dissipation under appropriate constraints. 

7. Appendix 

In this section, we shall derive the constitutive relations when the rate of dissipation is of the 
form 

£m(Dp( t) , D G ) = 7frD p (t).Dp(t) + ^ G D G .D G . (A.l) 

with the stored energy given by Eq. (3.16). Now, Eq. (3.15) reduces to 

(A.2) 
(A.3) 

(A.5) 



and 



T = pi + fjtBp(t) + %D P ((), 
T = AI + r ?G V K - 1 m D G V_ ) . 



Also, from Eqs. (A.2) and (A.3), 



(p-A) 



/' 



M B p(i))- 



Eqs. (|2~l5j), (l3~5j), we get 



Equating Eqs. (A.2) and (A.3), pre-multiplying and post-multiplying by V K , and using 

p(t) 

, v 
: V R m V 



(p-X)B pit) + fJ ,Bl (t) + r ^B 



»p(t) 



2 



pW 



(A.6) 



It follows that, the constitutive relations for the choices of Eq. (3.16) and Eq. (A.l) for the 
specific Helmholtz potential and the rate of dissipation, are 



T = P i + / iB p(4) + |v K -; (t) B p(t) v: 1 



p«' 



and 



Vg 



2 k pW B pW K p(t) 



I' 



tr(Bp(t))B p ( t ) + fiB 



2 

P(t) 



Vp-t 



B 



p{t)- 



(A.7) 
(A.8) 



Eq. (A.8) is the evolution equation of the natural configuration. 
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Figure 1 . Schematic of the natural configuration k v u\ corresponding to the cur- 
rent configuration K t and the relevant mappings from the tangent spaces of the 
same material point in kr, K t and K p (t) (above). The response from the natural 
configuration K p (t) is like a Kelvin- Voigt solid and the response of n p (t) from the 
reference configuration kr is purely dissipative. The corresponding one-dimesional 
spring dashpot analogy where a dashpot is in series with a Kelvin- Voigt element 
(below) . 
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FIGURE 2. Overall stretch of the current configuration from the reference config- 
uration (A) and square of the stretch from the natural configuration to the current 
configuration (B) as a function of non-dimensional time i for the creep experi- 
ment. For the loading process, T\\ = 1 and the unloading starts at t = 10. Plots 
for fj — 5, 10, 20 are shown. 
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FIGURE 3. Overall stretch of the current configuration from the reference config- 
uration (A) and square of the stretch from the natural configuration to the current 
configuration (B) as a function of non-dimensional time i for the creep experiment. 
For this case, the non-dimensional stress Tn = 5 and the ratio of the viscosities 
fj = 10. The unloading for the creep experiment starts at t — 10. 



16 




10 



15 



20 

t 



25 



30 



35 



40 



Figure 4. Overall stretch of the current configuration from the reference config- 
uration (A) and stretch from the natural configuration to the current configuration 
(\AB) as a function of non-dimensional time t for the creep experiment. For this 
case, the non-dimensional stress J\i = 1 and the ratio of the viscosities fj = 100000. 
Unloading for the creep experiment starts at t = 10. 
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FIGURE 5. Square of the stretch from the natural configuration to the current con- 
figuration (B), non-dimensional stress (Tii) plotted as functions of non-dimensional 
time i for various values of //, in the stress relaxation experiment. The initial con- 
dition for B was chosen as 2. 
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